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We study a classical ferromagnetic Heisenberg model in the presence of Dzyaloshinskii-Moriya 
interactions on the corner-shared triangle lattice formed by the Mn sites of MnSi. We show that a 
sizable spin helicity can be obtained only when the microscopic Moriya vectors lie parallel to the 
Mn-Mn bonds. Further, such vectors are shown to produce an unpinned helical order characterized 
by a particular ordering wavevector magnitude but unpinned direction, dubbed partial order, at 
physically realizable temperatures. A consequence of such an unpinned helical ordering is that the 
neutron scattering intensity is sharply peaked at this wavevector magnitude. The surface formed by 
connecting these wavevectors is a sphere, around which the neutron scattering weight is spread. We 
further show that the observed neutron scattering intensity can be anisotropic along this surface and 
that this anisotropy is dependent on the experimentalist's choice of lattice Bragg peak through a 
geometric factor. A neutron scattering measurement near the Bragg point (— , — ,0) naturally leads 
to a highest intensity along the (1,1,0) direction consistent with the observed anisotropy in MnSi 
[Pfleiderer et al. Nature 427 227 (2004)]. A possible mechanism for pinning the helical order, and 
a way to distinguish an ordered and a partially ordered state in the context of neutron scattering 
are discussed. 

PACS numbers: 75.10.Hk,75.25.+z,71.10.Hf 



Introduction: A recent paradigm for the discovery of 
new states of matter is the investigation of systems in 
close proximity to a quantum critical point (qcp). It 
was thus of great interest when high pressure neutron 
scattering^ of one of the most studied binary systems, the 
helimagnet MnSi, showed evidence of a new type of mag- 
netic state arising near a closely avoided qcp. Specifically, 
the periodicity of the spin helices was found to be fixed 
while their relative orientations had depinned, leading 
to a very disperse structure factor transverse to a well- 
defined radial ordering wavevector magnitude dubbed 
"partial order" . It has not been understood how such 
an anisotropic neutron scattering intensity arises from a 
partially ordered state. Starting from the same pressures, 
rather than exhibiting a quick recovery of Ap ~ T 2 as 
expected from spin fluctuation theory close to an itiner- 
ant qcp, resistivity measurements were found to exhibit 
a non-Fermi liquid character Ap ~ over an extended 
pressure range 2 ^. Moreover, the application of a small 
magnetic field, as might be expected to pin the helical 
axis, was found 2 -^ to lead to a low temperature recovery 
of Fermi liquid character, leading one to suspect an in- 
timate relation between the proposed partially ordered 
state and the unusual non-Fermi liquid resistivity preva- 
lent in zero field. 

It has long been believed that the essential physics 
present in itinerant hclimagnets on this lattice, in- 
cluding MnSi at ambient pressures, is the result of 
the interplay between itinerant ferromagnetic order and 
Dzyaloshinskii-Moriya (D-M) spin-orbit coupling terms. 
These tend to favor canted spin structures and are 



present because of the lack of inversion symmetry of 
the lattice^. A Landau-Ginzburg(LG) based phe- 
nomenological approach has led to many successes in- 
cluding the prediction of magnetic field induced wavevec- 
tor reorientations^. The LG free energy is expanded in 
terms of slowly- varying spin densities (S(r) in terms of 
which the DM interaction is given by S(r) • (V x S(r)). 
In such a continuum system, DM interaction has full ro- 
tational symmetry. As a result, the direction of ordering 
wave- vector is determined by a further anisotropic second 
order gradient term^ On the other hand, on a lattice, 
the DM interaction is given by Dy(S< x Sj), where Dy 
called the Moriya vector is defined on the link between 
two sites i and j. In general, the DM interaction breaks 
spin rotational symmetry due to a particular direction 
of Moriya vector, resulting in a preferential direction of 
the spin wave ordering wave- vector. The successful ap- 
plication of LG theory in the description of the physics of 
MnSi up to now lessened the motivation to study of the 
possible microscopic Moriya vectors of this lattice struc- 
ture. However, an understanding of the phase transition 
from an ordered helical state to the proposed partially 
ordered state as a function of applied pressure in this 
strongly correlated material requires a study beyond the 
standard LG theory^. 

More recently, to describe the partially ordered state, 
the addition of a pseudoscalar order parameter was 
argued^ to lead to a "blue quantum fog" with proper- 
ties of a chiral liquid. On the other hand, attempts to 
explain an observed shift of weight from the (111) di- 
rection of the ambient pressure ordered helimagnet, to 
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the apparent favoring of the (110) direction at high pres- 
sures, has led to theories which favor a weakly ordered 
state. This is argued to result from interactions between 
helimagnetic structures leading either to a new type of 
spin crysta&iS or a lattice of skyrmionsii. 

In this paper, we attempt to provide a microscopic 
model for helical ordering on the corner-shared triangle 
lattice formed by the Mn sites of MnSi, and to investigate 
one consequence of partial order in the context of the ob- 
served anisotropic neutron scattering intensity. We find 
that for reasonable spin-orbit couplings, the strongest 
component of the Moriya vectors must lie parallel to the 
Mn-Mn bonds, and favors nearly ideal simple helimag- 
netic structures with all spins perpendicular to the axis 
of propagation. However, in contrast to current folklore, 
at accessible temperatures the D-M interaction does not 
pin a direction of the helical ordering wavevector, but 
determines only its magnitude. A possible mechanism 
of pinning these helices to lie along the (1,1,1) direction 
at ambient pressure will be discussed. At high pressures, 
while the observation of high neutron scattering intensity 
along the (110) direction has stimulated ideas of crystal 
lattices featuring exotically ordered states, we show that 
such an anisotropic intensity can originate from a simple 
geometrical form factor, given a partially ordered state 
of unpinned helices. To confirm such an Occam's razor 
partial order further neutron scattering measurements 
relative to lattice Bragg peaks other than (— ,^-,0) are 
suggested^. 

Moriya vectors: The D-M interaction first came to 
wide attention to explain weak FM moments in largely 
antiFM a-Fe203 crystals by Dzyaloshinskiii^, but had 
been introduced several years earlier by Stevens^ where 
it was derived as a consequence of the inclusion of spin- 
orbit coupling in the Heisenberg model. Moriya^ in- 
cluded spin-orbit coupling in a consideration of the su- 
perexchange mechanism to provide the first microscopic 
derivation, introducing Dij (the Moriya vector) which 
is proportional to the first power of the spin-orbit cou- 
pling and an antisymmetric vector, and T ij — D '^ flli 

a symmetric tensor of second order in the spin-orbit cou- 
pling. For collinear spins, Yildirim et ali£ showed that 
the symmetric term can be of the same order of mag- 
nitude as the antisymmetric term, correcting Moriya's 
expression for the single U Hubbard model to read, 
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where [A <g> B\ llv = A^B^ and J,j is the regular Heisen- 
berg coupling. As the magnetic interactions in MnSi are 
thought to be predominantly ferromagnetic in origin, be- 
low we keep both terms linear and quadratic in , how- 
ever our main results do not appreciably change if we 
drop the contributions of the symmetric spin couplings. 
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FIG. 1: (Color online) Moriya's rules (a)-(e) for the deter- 
mination of the microscopic Moriya vectors^ in insulating 
materials due to superexchange interactions and (f) due to 
magnetic impurities at sites i and j in itinerant system a 21 ' 22 . 
Such vectors live on the bonds at the midpoint between atoms 
at sites i and j. Directions are assigned as shown when the 
crystal symmetries lead to: (a) the midpoint of Rij being a 
centre of inversion, (b) the midpoint of Rij being a mirror 
plane perpendicular to Rij, (c) a mirror plane parallel to Rij, 
(d) a 2-fold rotation axis perpendicular to Rij through the 
midpoint of Rij, (e) an n-fold axis parallel to Rij, and (f) 
an itinerant atom with a spin-orbit coupling living a distance 
Ri and Rj respectively from local moments at sites i and j. 
Symmetries (a) — (d) are not found for the MnSi lattice, which 
features a P2i3 symmetry. Case (e) applies, but only for ij 
being fourth nearest neighbours. 



To understand how to choose Moriya vector directions, 
it is useful to recall that MnSi is an itinerant magnetic 
system which shows features reminiscent of strongly cor- 
related systems. In particular, one finds an enhanced 
effective mass 10mJ^, a large effective moment at high 
temperatures (2.2 /is/Mni^, consistent with spin-| lo- 
cal moments at each site, above T c ), and a reasonable 
ordered magnetic moment ( 0.4 /Mn) 19 ' 20 . It has 
been argued that x-ray absorption measurements exhibits 
strong correlation effects^ Therefore, it is plausible that 
at low temperatures the electrons at each site have a 
dual character, producing both a small local moment 
and itinerant conduction electrons. In this case, magnetic 
RKKY interactions may arise between the local moments 
on each Mn site through these itinerant conduction elec- 
trons. The magnetic properties of such a model could find 
description in the language of the extended Heisenberg 
model. As the dominant magnetic interaction appears to 
be ferromagnetic in origin, as evidenced by the positive 
Curie offset of the magnetic susceptibility, we here con- 
sider the further simplification of keeping only the order 
1 terms of such an expansion, the nearest neighbour FM 
and DM interactions. One possible mechanism of helical 
ordering with wave-vector along (110) will be discussed 
below. 

In itinerant magnetic systems, it was shown that 
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FIG. 2: (Color online) a) The RKKY interaction of Fert and 
Levy chooses Moriya vectors parallel to the cross product be- 
tween the direction vectors from the local moment site to the 
spin-orbit site; b)-d) An orthogonal basis set of Moriya vec- 
tor choices consistent with the lattice symmetries: b) parallel 
to the bonds; c) parallel to the normal vectors, consistent 
with nearest neighbor RKKY contributions "RKKY I"; d) 
the cross-product between b) and c), consistent with further 
neighbor RKKY contributions "RKKY II" . 



the RKKY interaction is modified by the spin-orbit 
couplin g 21 ! 22 , generating microscopic Moriya vectors be- 
tween local moment sites. Their analysis was applied 
to non-magnetically doped CuMn x spin glass alloys in 
the dilute impurity limit. They found, that Z)y oc 

X r ! r 1 r^ i wnere ^ i s the spin-orbit coupling constant at a 
site of distance Ri and Rj respectively from two magnetic 
sites separated by distance Rij as shown in Fig. [T](f).For 
lattices like MnSi, where Moriya's rules as shown in Fig. 
[T] do not determine the direction of any nearest neigh- 
bor Moriya vectors, this allows one to construct a basis 
of directions consistent with the lattice symmetries, and 
potentially provides a mechanism for the realization of 
such an interaction. It is interesting to note that on the 
corner-shared tetrahedral lattice common to pyrochlorc 
systems the nearest neighbor RKKY interaction deter- 
mination agrees with the direction predicted by Moriya's 
rules 2 ^. 

Let us consider all possible sets of Moriya vectors con- 
sistent with the lattice symmetries. The lattice structure 
of the magnetic Mn sites is simple cubic with a four site 
basis, having an overall P2i3(T 4 ) symmetry. Sites can 
be labelled: {a, /3, 7, 6} — | — u, 1 — u), (1 — u, u+ 

5, 5 — u), (5— it, 1 — it, (u, u, u)} with UMn = 0.138. 

The silicon sites form a similar corner-shared triangle lat- 
tice structure with usi = 0.854. Fig. [2] (a) is a cartoon of 
the rule of Fig. Q] (f) as it would be applied to bonds be- 
tween local moments on neighboring Mn sites mediated 
by itinerant conduction electrons experiencing the spin- 



orbit scattering by a third Mn site. Basis vectors of Fig. 
[2tc) arise from considering the spin-orbit RKKY contri- 
bution from the third Mn site on each triangle, compo- 
nents along Fig. H{d) would arise from considering the 
closest Si atom as the spin-orbit site. An important lim- 
itation of our analysis is that the form of the Moriya 
vectors of Ref. ,21,23 was determined in the dilute impu- 
rity limit, whereas its application to MnSi would be in 
the dense impurity limit. As such, Fig. [2b) completes 
an orthogonal basis with Moriya vectors possibly of non- 
RKKY origin. While Moriya vectors as shown in Fig. [5] 
(b) are consistent with all Moriya rules and the lattice 
symmetries of MnSi, the microscopic derivation of such 
a possibility remains an interesting open problem. 

Partial order: For nearest neighbor ferromagnetic cou- 
plings (Jij < 0) between classical spins we have studied 
the low energy spectrum of the Hamiltonian given in Eq. 
[U We have used both mean fiel d 24 ' 25 and real space 
rigid rotor minimization 25 methods. The former method 
tends to over-estimate the likelihood of a spin structure 
to cant, but has no difficulty exploring extremely long pe- 
riod structures, and has been used to gain a qualitative 
feel for the relative strength of different direction cou- 
plings. In the latter method, we impose periodic bound- 
ary conditions on finite size structures with a given peri- 
odicity which allows us to access reasonably small order- 
ing wavevectors^ 5 -. 

In Fig. [3] (a), we show mean field results for the mag- 
nitude of the ordering wavevector for the three basis vec- 
tors. It is readily apparent that, given the size of the 
ordering wavevector of MnSi, the dominant contribution 
must arise from a choice of Moriya vectors parallel to the 
bond. This is confirmed by our rotor model calculations 
which find no evidence of canting when the Moriya vec- 
tors are chosen to lie in the primary RKKY I direction 
of Fig. [2] (c) , and internal structure and only very weak 
canting when they are chosen in the secondary RKKY II 
direction of Fig. f2] (d). In contrast, a choice of Moriya 
vectors parallel to the bonds, as shown in Fig. f2] (b), 
leads to almost perfect helimagnets with spin structures 
pinned in the plane orthogonal to the propagation direc- 
tion. 

Using a set of Moriya vectors lying parallel to the 
bonds, we investigate the directional energy dependence 
for helimagnetic wavevectors, Q . As shown in Fig. 
[3] (b) the spectrum along an infinite number of direc- 
tions, Q, is found to be degenerate within a precision of 
0(3 x 10 -6 J)2£. It is important to note that the wavevec- 
tors delineating members of this manifold lie on the sur- 
face of a sphere of radius Q | in momentum space. A 
pure helimagnet as, 

s(x) = sin(Q • ~~r)a + cos(Q • ~r^)S, (2) 

with ° x j?, = Q, where a and b are orthogonal unit vec- 

\axb\ 

tors is a reasonable approximation to the ground state 
spin structure as shown with open circles, lying only 
0(8 x 10~ 5 J) higher in energy. Mean field calcula- 
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FIG. 3: (Color online) (a) Ordering wavevector magnitude as 
a function of the D-M interaction for choices of the Moriya 
vectors shown in Fig. [21 where |Ad m | is set by the strength of 
\Dij\. Moriya vectors taken parallel to the bonds lead to sim- 
ple spirals which can account for a period as short as MnSi; 
those of RKKY I do not cant spins; while those of RKKY 
II lead to internal structure and very weak canting, (b) 
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Q = - v/gjr+9p + gf an d Moriya vectors parallel to the 
bonds, the dispersion is to a very good approximation given 
by a degenerate sphere in momentum space. Filled symbols: 
rigid rotor solutions, solid lines along qOO, qqO, qqq: mean 
field, and open circles: pure helimagnet. (c) At low temper- 
atures radial cuts of the structure factor can become sharply 
peaked along any direction in comparison to transverse cuts 
for energetic reasons as discussed in the text. 



tions (which feature a soft spin constraint) are found 
to provide a reasonable approximation to the energet- 
ics of the true minimized spin structures at small Q 
as shown by the solid lines of Fig. [3] (b), allowing 
a smooth interpolation of the low energy states avail- 
able to a finite temperature system. This means that 
we can approximate the probability of a given state as 



probability, as seen in Fig. [3J (c), is the dominant con- 
tributor to the radial structure factor. 

Anisotropic structure factor: An interesting question 
is whether the anisotropic neutron scattering intensity 
strongly favoring the (110) direction under pressure in 
MnSi can be understood to arise from a partially or- 
dered state, that is, from an unpinned helical order. 
Here we show that an anisotropic intensity can be un- 
derstood to originate from a geometric factor present 
due to the dipolar nature of the neutron-electron in- 
teraction. The interaction between the magnetic mo- 
ments of the neutrons and electron spin leads to the form 
Hi n t(sj) oc e~ tK ' Rj jt n ■ (Q x ~~Sj x Q) where ~ft n is the 

neutron magnetic moment, HQ the momentum transfer, 
and s j the spin at site j. The elastic contribution to the 
differential cross-section is found in the Born approxima- 
tion to be^l, 
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where a, (3 run through {x, y,z}. In Eq. [2 we found 
that a reasonable approximation to the ground state was 
a pure spin spiral with a small wavevector, Q. As the 
measured MnSi datai lie about the (110) Bragg peak we 

have Q — Q + (^, ^,0), and one expects that both 
an atomic form factor and an angular dependence should 
arise. Assuming the form of Eq. [3 at zero temperature, 
we find the structure factor to be proportional to, 



(l-cos(87ru))(l-i((a-Q) 2 



(4) 



Note that the geometric form factor in Eq. [3] indicates 
that unpolarized neutron scattering is insensitive to the 
component of a spin structure parallel to the momentum 
transfer of the scattering neutrons. Hence, for a truly 
degenerate simple helimagnet measured with respect to 
the (0,0,0) lattice Bragg peak one expects to see no an- 
gular structure, while measurements with respect to the 
(1,1,0) lattice Bragg peak will generate diffuse intensities 
over the surface of a sphere in momentum space which 
peaks along (1, 1,0). This is illustrated in Fig. [U 

Discussion: The nature of the high pressure state 
of MnSi remains a topic of active experimental debate. 
While some have argued that MnSi undergoes a first- 
order phase transition at quite low pressures^, the pre- 
dominant focus has been on understanding the detailed 
magnetic properties of the observed signature of partial 
order. While nuclear magnetic resonance studies of pow- 
dered samples found decreasing magnetic moments well 
above the critical pressure above which the helices have 
dcpinned3^, a careful muon spin relaxation study^ of 
single crystals did not find ordered moments. Rather, 
Uemura et al. 30 argued for a "slow but dynamic order" 
as responsible for the diffuse neutron scattering signal 
observed, suggesting that the spins must fluctuate on a 
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FIG. 4: (Color online) (a) While the spin plane (aqua) of the 
helimagnet pins perpendicular to the ordering wavevector Q, 
this spin plane may not be perpendicular to the momentum 
transfer of the neutrons, Q = Q + (^-,^ L ,0). Note that con- 
structive interference leads to a magnetic signal at all points 
related to Q by a reciprocal lattice vector, (c)-(d) Variation of 
the structure factor arising due to the geometric factor (Eq. Q 
about the lattice Bragg peak (— , — ,0) assuming a pure spiral 
ground state of the form Eq. [2] with |Q|=0.043A. 8 and <j> are 
chosen such that Q — |Q|(sin((9) cos(<j!>),sin(#) sin(0),cos(f?)). 
Note the maximum occurs at (1,1,0) and (-1,-1,0) but not at 
(1,0,1), (0,1,1). 



a sizable contraction of the lattice has been observed by 
neutron Larmor diffractior*2I to be concomitant with the 
onset of non-Fermi liquid behaviour at high pressures, 
ruling out most quantum critical scenarios (ex. magnetic 
rotons 3 ™) . It would appear that this lattice contraction is 
responsible for the depinning of the helimagnetic order. 



In this light, we briefly comment on the pinning of 
the ambient pressure (1,1,1) ordered state. The local 
environment of the Mn site is surrounded by 7 Si atoms, 
with the shortest Mn-Si distance along the local (1,1,1) 
axis. It therefore seems plausible that crystal/ligand field 
effects might produce a local (1,1,1) anisotropy of the 
spin structure. One term which might arise from a full 
consideration of spin-orbit interactions in this context is 
of the form, — (Yli v( s iising) 4 wnere the sum runs over 
the number of unit cells and the four sites within the 
unit cell, ~s" Ising — ( s i • ~e v )e v and ~e ' v is the unit 
vector associated with the local (1,1,1) anisotropy at site 
v. Assuming this term is small, we can evaluate the 
energy change such a term would introduce to the simple 
helimagnet. For £ > this term is found to lower the 
energy of the 100, 110 and 111 directions by -6£, -6.5£ and 
-8C respectively per unit cell, effectively favoring (1,1,1). 
The increase of pressure would likely increase the overlap 
of the various Si atoms with the Mn sites, reducing the 
strength of pinning. 

We have shown that the geometric form factor provides 
one possible explanation of the higher intensity along the 
(1,1,0) direction measured by neutron scattering 1 about 
the (— , — ,0) Bragg peak. However, we cannot rule out 
the possibility that additional terms may be present and 
weakly favor ordered helices along (1,1,0). In particu- 
lar, we have showr*22 that the addition of antiFM 2nd 
and 3rd nearest neighbors generically leads to a (1,1,0) 
ordered helical state. We expect that this (1,1,0) or- 
dered helical state will be robust in the presence of the 
DM interaction, as the DM interaction does not pin the 
direction^ To distinguish these two different proposals, 
we propose further neutron scattering around different 
Bragg points such as (2^,0,0) is desirable-in the par- 
tially ordered state a higher intensity moves to (1,0,0) 
while in a (1,1,0) ordered helical state, the neutron scat- 
tering intensity remains peaked at (1,1,0). Even within 
the existing study, it is our belief that the "rotary electric 
shaver" picture of the intensity created by mathematica 
for publication by Pfleiderer et alX, while consistent with 
their data, may have inadvertantly symmetrized their in- 
tensities, as we would not expect a similar enhancement 
near (0,1,1) and (1,0,1) arising solely from the partial 
order here considered. 

Summary: To conclude, we have shown that the domi- 
nant contribution of the Moriya vectors in MnSi must lie 
parallel to the Mn-Mn bonds. This leads to a dispersion 
minimum at a particular helix pitch independent of the 
direction of order. The corresponding momentum space 
sphere has a much smaller energetic spread than acces- 
sible temperatures, leading to a sharply peaking radial 
structure factor, and a very disperse transverse structure 
factor. We offer a natural explanation for the higher 
intensity of the structure factor at (1,1,0) as observed 
experimentally. The application of weak magnetic holds 
would be expected to break this degeneracy, pin the ori- 
entation of the helimagnets, and it is tempting to argue 
that the pinning of the helical order direction leads to 



6 



the recovery of Ap ~ T 2 behavior seen. 
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[Note Added] After a completion of our work, the au- 
thor of Ref. HI agreed that that magnetic scattering near 
a nuclear Bragg spot is subject to a selection rule, which 
has not been discussed in Ref. 01 However, he pointed 
out that Fig. 2 (b) and (d) in Ref. [l| show that the inten- 
sity drops off much faster than what one would expect 



from the effect of geometrical form factor in unpinned 
helicial state. While we show that the intensity drops 
to 50%, when the azimutal angle, <j> becomes ±7r/2 at a 
fixed 9 = ir/2 as shown by a line with circles in Fig. 0](b), 
the experimentally observed intensity drops to almost 
near ±7r/6, indicating neither an unpinned helix nor 
an ordered state. Other (110) lattice points exhibit the 
maxium intensity which have not been reported in Ref. 
1. These confirm that further studies such as neutron 
scattering around other Bragg lattice points are needed 
to determine the detailed nature of the partially ordered 
state. 
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